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Abstract. The method originally developed for the exact calculations in QED theory 
is applied for the calculation of NLO effects in QCD Compton processes. QCD correc- 
tions to the structure functions and sum rules are obtained. Different interpretations 
of the NLO effects due to finite quark mass are discussed. 
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1 Introduction 

A lot of articles dedicated to NLO QCD corrections to polarized DIS have been 
published. Calculations in most of them are performed in an assumption that 
the quark mass is equal to zero (see [1, 2, 3, 4, 5, 6, 7] and review [8]). However 
there are some works when authors estimate a finite quark mass effect for both 
unpolarized [9] and polarized structure functions [10, 11, 12, 13, 14, 15, 16]. The 
main difference of these two approaches is the method of tending the quark mass 
to zero and the procedure of integration of the squared matrix element over the 
phase space of emitted gluon. In the massless approsch a fermion mass is equal 
to zero before the integration while in the massive one this quantity goes to zero 
after the integration and survives only in LO terms. 

Notice that QCD and QED radiative corrections (RC) having different origins 
possess some common features on the one-loop level. If our consideration is 
restricted to so-called QCD Compton process then both corrections should be 
described by the identical sets of Feynman graphs. The transition from the strong 
radiative effects to the electromagnetic ones could be performed by the following 
replacement: 

-a s -> e 2 q a QE D- (1) 

In the present paper we investigate QCD Compton process for polarized 
nucleon target by the method developed for exact calculations [17, 18, 19, 20, 
21, 22] in QED theory for massive polarized fermions but never used for QCD. 
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We estimate the influence of the finite quark mass effects on the value of the 
polarized sum rules and discuss some different interpretations of the first moment 
of the polarized structure function gi . 

The present paper is organized in the following way. The method of calcula- 
tion is described in Section 2. The explicite expressions for the QCD-improved 
structure functions as well as their first moments can be found in Section 3. 
Then some conclusions are made in Section 4. We also discuss there different 
points of view how the NLO finite quark mass effects influence the first moment 
of f/i. Details of the integration procedure are presented in Appendix. 



2 Method of Calculation. 

The cross section of polarized lepton-nucleon DIS 

l(kuZ) + N(p,rj)^l(k2)+X (2) 

on the level of the one-photon exchange can be presented as a convolution of the 
well-known lepton tensor: 

= ^Sp 7 M (fci + m)(l + 75£)7^(fo + m), (3) 

and the hadronic one: 

W IW = -g^F 1 + ^F 2 
pq 

+—e^q X (v' J 9i + {rf - ^V).9 2 ). (4) 
pq * pq 

Here p and r\ are an initial 4-momentum of the target and its polarization vector 
respectively and q is a 4-momentum of the virtual photon (p 2 = M 2 , q 2 = —Q 2 ). 

On the Born level in the frame of the naive parton model the structure 
functions contributed to the hadronic tensor W® v have the following form: 

9 

F°(x) = xJ2e 2 q f q (x), 
i 

9 

g° 2 (x) - 0. (5) 



Here f g (x) and Af q (x) are unpolarized and polarized parton distributions re- 
spectively 1 , and x = Q 2 /2pq is the standard Bjorken variable. 

1 Hcrc and later the independence of the polarized parton distribution on the polarization 
vector is assumed. See discussion in [5]. 



http : //link. springer .de/link/ service/ journals/ 10105/ index.html 



EPJdirect C5, 1-12 (2008) Springer- Verlag 



3 



The one-loop lowest-order RC to the hadronic current consists of two parts 
with contributions calculated in different ways: 

W£ 1 «" = WV+W*. (6) 

The first part of the one-loop RC to the hadronic current comes from the gluon 
exchange and it will be described later. The second one corresponds to the gluon 
emission and requires the integration over its phase space: 



W, 



x[v$?f q (x/z) + w$Af q (x/z)], (7) 



where 



= Sp r^aiplg +m q )r av {p2q +' 



w 



= Sp r^ a {p lq +m q )i 5 fir a „(p 2q + m q ), (8) 



and p\ q (p2q) is an initial (final) 4- momentum of the quark {p\ q = p\ q = m q ) 
Quantities r^ a and T av are defined as 

^)tQ — ^^q 







Vq 


Uq 






Vq 


Uq 



r av = 2(1^-^^-^^, (9) 

where 

Q Plq _ Pjq^ (1Q) 

Vq U q 

and v q = 2kp\ qi u q = 2kp 2q - A variable z = Q 2 /(Q 2 + u q ) is a standard one used 
for the description of QCD effects. 

Within the requirements of the naive parton model in (7) we assume that 
the quark mass depends on the variable z as 

m q = x/zM (11) 

for bremsstrahlung process and 

m = m q (z = 1) = xM (12) 

for non-radiative one. Notice that the factor before M in these formulae coincides 
with the argument of the corresponding parton distribution. 

Both contributions in right hand side in (6) include the infrared divergences 
to be carefully considered for canceling. Similar to QED we use an identity: 



W = W H - W 1H + W 1H = W b + W 1H (13) 
Here an infrared part can be written in the following way: 

K? = t^ SIRw °^ ( 14 ) 
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where is a contribution of g-quark to the hadronic tensor on the Born level 



and S IR should be decomposed in a sum of a soft and hard parts: 
s ir = if <Pk F iR = I f ^ F m e{e _ ko) 

7T J k a IT J Kq 

+ 1 J ^ FlRe ^ £ ) = % + (15) 

Here ko is the energy of the emitted gluon, e is an infinitesimal parameter of the 
separation and 

F IR = + ^ (91+^1 - H$ _ -1 1 ,;,.(,) 



(Q 2 + Uq - V q ) \ UqVq U 



Notice that there is some arbitrariness in the definition of F IR : the asymp- 
totic expression (k — > 0) is fixed. We choose it like (16) to write the QCD- 
improved structure functions (see formulae (27)) in the simplest form. Recall 
that the quark mass in (16) is defined by (11). 

The way of calculation of typical integrals in the soft and hard parts of 8 IR 
is widely discussed in [21, 22]. After explicit integration the sum of the soft 

5 IR ft = 2(l 1) (V* + log ^ - ll + l + 1 - tt 2 /6 (17) 

and the hard parts 

S IR rd = 2(l 1) log ^ + hi + l l v - \ll -l -l v - tt 2 /6 (18) 

gives the result independent on the parameter e. Here 

Q 2 1 - x 

Z = log— o, Z„=log . (19) 

mg x 

The pole term which corresponds to the infrared divergence is contained in 

p "' = ^- i + b° +l °^ < 20 » 

The arbitrary parameter /n has a dimension of a mass and 7_e is the Eulcr con- 
stant. 

To extract some information about QCD contributions to the polarized struc- 
ture functions gi and gi an integration in W^ v over the gluon phase space should 
be performed without any assumptions about the polarization vector r\. So the 
technique of tensor integration has to be applied in this case. All exact expres- 
sions for tensor, vector and scalar integrations can be found in the appendix. 
The ultrarelativistic limit m q — > can be performed in a straightforward way, 
with the exception of two terms: f(x)/r, u q f(x/z)/r 2 . The first one including 
1/r appears only in the substracting parts of W^ v . Since the corresponding par- 
ton distributions do not depend on u q it can be integrated over u q analytically. 
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The second one u q /r 2 comes from W^. The ultrarelativistic limit in can 
be applied after the transformation 

J du q ^f(x/z) = J du q f{x/z) u ~ f{x) + f{x) J du q ^ (21) 

and the analytical integration of the last term. Here / is any (either unpolarized 
or polarized) parton distribution. The first term in the right hand side of (21) 
gives the contribution to W^ v . The second one together with the result coming 
from the integration 1/t over u q gives an additional factorized correction S m : 

S m = -I(3i + 3/„ + l). (22) 
An explicit expression for the gluon exchange contributions reads: 



where 



3 7T 



, : _ 2 (p^ + log^)(Z -l)-^ + I/ - 2 + - (24) 



is a vertex correction part factorizing before W®%, and 

iM 
pq 



K» MM = & E ^ W [rf - f/] E ASM (25) 



is the quark anomalous magnetic moment. 

The infrared terms P IR , parameters /i, e and the squared logarithms con- 
taining mass singularities are canceled in the sum 

g IH + S rn + § V = 1 L lQlv _ £ + _ 7^ _ 5 _ ^ 



= o<5 9 (26) 



that is identical to (C.17) from [23]. 

Since the result of the analytical integration has the same tensor structure 
as the usual hadronic tensor in polarized DIS, the coefficients before the cor- 
responding tensor structures (like g^, Pp.p v ... ) in the sum W^ v + W^ u can 
be interpreted as one-loop QCD contributions to the corresponding structure 
functions. 



3 QCD-improved structure functions 

In the last section we obtained one-loop correction to the hadronic tensor and 
structure functions. The explicit result for them is: 

Fi(x,Q 2 ) = l.[F 2 (x,Q 2 )-F L (x,Q% 
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F2(x,Q 2 ) 

Fl(x,Q 2 ) 
9i(x,Q 2 ) 
92(x,Q 2 ) 



4a. 
"3^ 



1 

J24 dzf q (x/z), 
^e 2 q Af q (x,Q% 

^£e*{(l-2Z ? -log(l-aO)4f g (*) 



(4l g - 4\ogz{l - z)- 12 



Afqjx) 



-)Af q (x/z) + 



} 



(27) 



The Q -dependent unpolarized and polarized parton distributions are denned as 

l 







x (i, - log z(l - z)) - + 3z + 4) /, (H) 



Af 9 (*,Q 2 ) = (1 + |^)A/ 9 M + / y 



.1 + z 2 
1 - z 



x (i, - log z(l - z)) - - — + Az + l) Af q (-J 



(28) 



where l q = \ogQ 2 /m 2 and S q is defined by (26). 

Notice that our formulae (27,28) are in an agreement with ones obtained 
earlier. The unpolarized structure functions coincide with (2.24,2.49) from [9], 
gi(x, Q 2 ) corresponds to the expression (13) in [10]. At last g 2 {x,Q 2 ) can be 
compared with the sum (18) and (19) from [11]. 

Let us calculate the difference between QCD-improved parton distributions 
found in both of the discussed schemes: 

|^(/r =0 (*> q 2 ) /r^v, q 2 )) = [{p ir 

+ log^)^ + 2l v )-l 2 v -l v + l\f q (x) 
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+ 



/f[(V 



: ^^ + i og !^ + i og (i 



-OP' 



+ 



1 



z 



—z 
2 



yU Z 



+ 



3tt 

4a., 



(Z\/!"-°(x,Q 2 ) - Af™«*°(x,Q 2 )) = \(P 



jIR 



too , .3 „, „ 
+ 1 ° g M K 2 + " } 



ly ly ~\~ 1 



/^rfi±f! ( p« + i og !^ 



4f«(z) 
log(l - z)) 



+ 



1 3 

l-z ~ 2 Z ~ 



-^-(p^ + log^+logi 

- Z V /Lt 



+ 



\)Af q {x) 



(29) 



It can be seen from (29) that the partonic distributions within the two schemes 
have the same structure and all leading logarithmic terms are canceled after the 
following replacement: 

logm,, -> log /i-P Jfl . (30) 

The standard QCD results include infinite terms P IR (see ref.[6]) that corre- 
sponds to quark mass singularity terms in our formulae. These divergent and 
scheme dependent terms in the classical results arc rejected by the additional 
procedure of the renormalization of the parton distributions (see ref.[8] for de- 
tails). Contrary, our formulae (28) being dependent on logm 9j o do not require 
any renormalization. The main origin of differences of the non-leading logarith- 
mic terms in these formulae comes from the calculation procedure itself since 
the integration over the emission gluon phase space and tending of the quark 
mass to zero are not commutative with each other. 

In order to find one-loop QCD contributions to the sum rules it is necessary 
to integrate (27) over the scaling variable x. After an identical recombination 
like 



dx J dzR(z)f q (x/z) = 

X 

an explicit integration over z gives 



dzzR(z) J dxf q {x) 
o 



(31) 



1 1 

J dxF^x.Q 2 ) = (l - ^) J dxF°{x) 



dx 



F 2 (x,Q 2 ) 



dx 
x 
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l l 






1 



dxg 2 (x,Q 2 )= 0, 



(32) 



o 



where the structure functions with index "0" are defined by (5). 

Notice that QCD RC to the first moment to the unpolarized structure func- 
tions as well as gi have identical values for both massive and massless quark 
approaches. At the same time these two methods give the different results for 
the QCD corrections to the first moment of g\. 

Till present time the calculations have been performed within the naive par- 
ton model when the quark mass is defined as in (11, 12). Another possibility 
is to consider it as a constant like it was assumed in [24]. An additional factor 
can depend only on z and x. So the result for the structure functions would 
differ only in an argument of logarithm containing the quark mass. An addi- 
tional contribution to the sum rules can appear after integration over x and z. 
However it could be shown by the explicit calculations that this contribution is 
equal to zero. Indeed an additional integrand has DGLAP-like structure, so this 
cancellation is consequence of DGLAP equations. Thus the final results for sum 
rules in these two cases are identical. 

4 Discussion and Conclusion 

In this paper we applied the approach traditionally used in QED and electroweak 
theory for calculation of QCD correction to the DIS structure functions and sum 
rules. Since the quark was considered massive within this approach, it allowed 
us to estimate the finite quark mass effects at the NLO level. LO correction 
contributes to the DIS structure functions but vanishes for the sum rules, so this 
mass effects are important just for the sum rules. 

We found that there is no any additional effect for the first moment of the 
unpolarized structure functions as well as for g 2 . However there is some non-zero 
correction to the first moment of g\ and as a result to the Ellis- Jaffc and Bjorkcn 
sum rules. A value of the correction is in agreement with the results of refs.[10, 
13], obtained by different methods. We confirm also the statement of [10] that 
the classical value of correction to the first moment of g\ is reproduced if we take 
into account the leading term of the polarization vector r\ — p\ q /m q . However 
contrary to the paper we would interpret the result with (-5/3) in cqs.(32) as 
physical one. In our calculation we took the proton (quark) polarization vector 
in a general form. Using of the exact representation of the vector [18] (k\ is an 
initial lepton momentum and k\ = mf): 



leads to exactly the same result (with -5/3) for the first moment of g\. The 
second term giving the additional contribution (-2/3) cannot be neglected within 




1 



(33) 
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the approximation under consideration. It can be easyly understood from pure 
calculation of QED and electroweak corrections to the lepton current in DIS 
process [20, 25], where the lepton polarization vector looks like (33), and this 
additional contribution was extracted to separate non-vanishing contribution. 
The result with (-5/3) was obtained under similar assumptions in the report 
[12] . The additional correction was discussed both within OPE (operator product 
expansion) and within improved quark-parton model in paper [13]. It was shown 
in this paper that there is no contradiction between this result and classical 
correction obtained for massless QCD, if we carefully take into account the finite 
quark mass effects for coefficient function and matrix clement. As it was reviewed 
in recent paper [15] the renormalization of the axial-vector current completely 
suppresses the additional contribution obtained within massive approach (see 
section 4 of the paper). 

We note that Burkhard-Cottingham sum rule is held withing taking into ac- 
count mass effects at the NLO level. The target mass correction of the next order 
(~ uiq/Q 2 ) was analysed and was found negative in the paper [14]. However, in 
the case of conserved currents (see [16]) the polarized structure function g 2 with 
the target mass correction obeys not only Burkhard-Cottingham sum rule but 
Wandzura-Wilczek relation too. 

Estimating the diagrams with one gluon radiation we have also result for 
QED radiative correction to hadronic current. The calculation within the quark 
parton model shows that QED correction is not so large, however there are at 
least two arguments to be taken into account. First, the DIS structure functions 
are defined in the one-photon exchange approximation as an objects including 
only strong interaction. It means that transferring from the cross section to the 
structure functions we neglect these QED effects, so their contributions have to 
be included to systematic error of corresponding measurements. Second reason 
is that there exist some measurements where QED correction is important. One 
of the examples is calculation of a s by comparing theoretical and experimental 
values for the Bjorken sum rule [26]. In this case a s is suggested to be extracted 
from QCD corrections to the Bjorken sum rule, which is a series over a s , and at 
least five terms of the expansion are known. Simple estimation shows that QED 
contribution is comparable with fourth (or even third) term of the expansion. 

There are several papers devoted to calculation of QED and electroweak 
corrections to the hadronic current [24, 25, 27]. As usual methods similar to 
ours are used for that. Positive moment here is the keeping quark mass non- 
zero. From the other side there are some effects which are not considered within 
QED calculations: taking into account the confiment effect in integration of soft 
region, consideration of the quarks as non-free particles. Thus the methods (see 
for example [12]) developed for careful treatment of the QCD effects should be 
applied for the analysis of photon emission from the hadronic current. However, 
probably the best way to take into account the photon radiative correction to 
the hadronic current is further generalization of OPE technique to include QED 
effects. 
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Appendix 



The gluon emission phase space on the parton level can be written in a covariant 
form: 

d 3 k {Q 2 + U q — V q ) d,Uqdv q dzi 

k ' 2{Q 2 + u q ) y /-A(k 1 ,k 2 ,pi q ,P2 q )' 
Here v q = 2kp\ q , u q = 2kp 2q , Z\ = 2kki, A(k\, k 2 ,Pi q ,P2 q ) is Gram determinant, 



(A.l) 



(Q 2 



ArriqQ 1 and p\ q (p 2q ) is an initial (final) quark momentum. 



Since the gluon is radiated from the quark legs it is clear that the matrix 
element squared for this process does not depend on z\. Therefore an analytical 
integration of the right part of (A.l) over z\ can be performed 



/ 



dz\ 



2 v /-A(k 1 ,k 2 ,pi q ,P2 q ) 



(A.2) 



where the limits z ™ ax l min are defined as the solutions of the equation 

A{k\, k 2 ,Pl q ,P2q) = 0. 

The following expression 



q 

J[A] = J dv q i 



(A.3) 



(A.4) 



for the integration over v q will be used in this appendix. Explicit expressions for 
the limits of integration which can be obtained as the solution of the equation 



~?nax ,mm 

Zi — Z, 



read: 



max/min 



U q (Q 2 +U q + 2m q ± y/\q) 



2t 



(A.5) 



i q , ,.„ q . 

Tensor and vector integrals can be reduced to scalar ones using the following 



expressions: 



kii k u A 



J 



k^A 



2\„ 



T\A 



X 2 p^Pv 



T 2 A 



V 2 A 



T 3 A 
V 2 A 



T 4 A 



(A.6) 
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The coefficients 7* and Vi are defined as the solution of the systems of equations 
obtained by the convolutions (A. 6) with p\ q , q and : 

Ti = (u q - v q ) 2 m 2 + u q v q (v q - u q - Q 2 ), 

T 2 = ((U q + Vg) 2 + 2u q V q )Q 4 + 2( Uq + 2Vq)Q 2 U q 

+ (u 2 q -2m 2 q Q 2 )(u q -v q ) 2 , 
T 3 = {Au q -v q )Q 2 u q v q + {2u q + v q )Q 3 v q 

-m 2 (u q - V q )(3u q (Q 2 +U q - V q ) + Q 2 V q ), 

T 4 = (Q 2 + u q ) 2 v 2 q +6{u q -v q ) 2 m q 

-2m 2 v q [(Q 2 + u q )(3u q - 2v q ) - u q v q ], 
Vi= Q 2 (u q + v q )+u q (u q -v q ), 

V 2 = v q (u q + Q 2 ) + 2m 2 q (v q -u q ). (A.7) 
The list of scalar integrals over variable v q in NLO approximation reads as 



V 

q 


m 2 u q ' 


Vq 


(Q 2 +u q + 2m\ + ^f\) 2 
4m q T 


J[l] = 


r 


J[v g ] - 


u 2 q (Q 2 + u q + 2m 2 q ) v % 
2t 2 



(A.8) 



References 

1. G. Altarelli, G. Parisi: Nucl. Phys. B126 (1977) 298 

2. J. Kodaira, S.Matsuda, K. Sasaki T.Uematsu: Nucl. Phys. B159 (1979) 99 

3. J. Kodaira: Nucl. Phys. B165 (1979) 129 

4. J. Kodaira, S. Matsuda, M. Muta, T. Uematsu, K. Sasaki: Phys. Rev. D20 (1979) 
627 

5. M. Anselmino, A. Efremov, E.Leader: Phys.Rept. 261 (1995) 1 

6. P. Ratcliffe: Nucl.Phys. B223 (1983) 45 

7. M. Gliick, E. Reya, M. Stratmann, W. Vogelsang: Phys. Rev. D 53 (1996) 4775 

8. B. Lampe, E. Reya: Phys.Rept. 332 (2000) 1 

9. B. Humpert, W.L. van Neerven: Nucl. Phys. B184 (1981) 498 

10. R. Mertig, W.L. van Neerven: Z. Phys. C C60 (1993) 489: Erratum-ibid. C 65 
(1995) 360 

11. G. Altarelli, B. Lampe, P. Nason, G. Ridolfi. Phys. Lett. B334, 184 (1994). 

12. O. Teryaev: Contribution to International Symp. Dubna Deutreron '95, Dubna 
Russia 1995, hep-ph/9508374. 

13. O. Teryaev, O. Veretin: hep-ph/9602362 

14. I.V. Musatov, O.V. Teryaev, A. Schafer: Phys.Rev. D 57 (1998) 7041 

15. V. Ravindran, W.L. van Neerven: hep-ph/0102280 

16. J. Bliimlein, A. Tkabladze: Nucl. Phys. B 553 (1999) 453 



http : //link. springer .de/link/ service/ journals/10105/ index.html 



EPJdirect C5, 1-12 (2008) Springer- Verlag 



12 



17. D. Bardin, N. Shumeiko: Nucl. Phys. B 127 (1977) 242 

18. I. Akushevich, N. Shumeiko: J. Phys. G 20 (1994) 513 

19. I. Akushevich, A. Ilyichev, N. Shumeiko, A. Soroko, A. Tolkachev: Comp. Phys. 
Commun. 104 (1997) 201 

20. I. Akushevich, A. Ilyichev, N. Shumeiko: J. Phys. G 24 (1998) 

21. A. Akhundov, D. Bardin, L. Kalinovskaya, T. Riemann: Fortsch. Phys. 44 (1996) 
373 

22. D. Bardin, L. Kalinovskaya: DESY-97-230 Dec 1997. 54pp., hep-ph/9712310. 

23. D. Bardin, C. Burdik, P. Christova , T. Riemann: Z. Phys. 42 (1989) 679 

24. I.Akushevich, A. Ilyichev, N. Shumeiko: Phys. At. Nucl. 58 (1995) 1919 

25. D. Bardin, J. Bliimlein, P. Christova, L. Kalinovskaya: Nucl.Phys. B 506 (1997) 
295 

26. J. Ellis, M. Karliner: Phys. Lett. B 341 (1995) 397 

27. H. Spiesberger: Phys. Rev. D D52 (1995) 4936 



http : //link. springer .de/link/ service/ journals/10105/ index.html 



